Properties of tunnel Josephson junctions with a ferromagnetic interlay er 
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We investigate superconductor/insulator/ferromagnet/superconductor (SIFS) tunnel Josephson junctions in 
the dirty limit, using the quasiclassical theory. We formulate a quantitative model describing the oscillations of 
critical current as a function of thickness of the ferromagnetic layer and use this model to fit recent experimental 
data. We also calculate quantitatively the density of states (DOS) in this type of junctions and compare DOS 
oscillations with those of the critical current. 
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I. INTRODUCTION 

It is well known that superconductivity and ferromagnetism 
are two competing orders, however their interplay can be real- 
ized when the two interactions are spatially separated. In this 
case the coexistence of the two orderings is due to the prox- 
imity effect''^-^. Experimentally this situation can be realized 
in superconductor/ferromagnet (S/F) hybrid structures. The 
main manifestation of the proximity effect in S/F structures is 
the damped oscillatory behavior of the superconducting cor- 
relations in the F layers. Two characteristic lengths of the 
decay and oscillations are, correspondingly, t,f\ and E,f2. Un- 
usual proximity effect in S/F layered structures leads to a num- 
ber of striking phenomena like nonmonotonic dependence of 
their critical temperature and oscillations of critical current in 
S/F/S Josephson junctions upon the F layer thickness. Neg- 
ative sign of the critical current corresponds to the so-called 
n state. Spontaneous n phase shifts in S/F/S junctions were 
observed experimentally . ■*'^A7,8,9,10,11,12,13,14,15 

SIFS junctions, i.e. S/F/S trilayers with one transparent in- 
terface and one tunnel barrier between S and F layers, repre- 
sent practically interesting case of tt junctions. SIFS structure 
offers the freedom to tune the critical current density over a 
wide range and at the same time to realize high values of a 
product of the junction critical current Ic and its normal state 
resistance R^}'^'^^ In addition, Nb based tunnel junctions are 
usually underdamped, which is desired for many applications. 
SIFS K junctions have been proposed as potential logic ele- 
ments in superconducting logic circuits."' SIFS junctions are 
also interesting from the fundamental point of view since they 
provide a convenient model system for a comparative study 
between O-tt transitions observed from the critical current and 
from the density of states (DOS). At the same time, despite 
such an interest, there is no complete theory yet of SIFS junc- 
tions which could provide quantitative predictions for critical 
current and DOS in such structures. All existing theories dealt 
only with a number of limiting cases, when either linearized 
quasiclassical equations can be used for analysis'^ (e.g. tem- 
perature range near critical temperature, small transparency of 
interfaces) or thickness of the F layer is smalli^ compared 



to the decay characteristic length ^fi. Further, in symmetric 
S/F/S junctions, the extension of theory to the case of nonho- 
mogeneous magnetization and large mean free path was per- 
formed in Refs. 3. 18. 

The purpose of this work is to provide a quantitative model 
describing the behavior of critical current and DOS in SIFS 
junctions as a function of parameters characterizing material 
properties of the S, F layers and the S/F interface transparency. 
The model provides a tool to fit experimental data in existing 
SIFS junctions. 

The paper is organized as follows. In the next section 
we formulate the theoretical model and basic equations. In 
Sec.|lII]we solve nonlinear Usadel equations, apply solutions 
for calculation of critical current in SIFS junctions with long 
ferromagnetic layer, df ^ E,fi, and fit recent experimental 
data. In Sec.|IV]we perform numerical calculations for critical 
current in a SIFS junction with arbitrary length of the F layer. 
In Sec. [V] we numerically calculate DOS in the ferromagnetic 
interlayer, and then summarize results in Sec. |VT] 



II. MODEL AND BASIC EQUATIONS 

The model of an S/F/S junction we are going to study is de- 
picted in Fig. [Hand consists of a ferromagnetic layer of thick- 
ness df and two thick superconducting electrodes along the 
X direction. Left and right superconductor/ferromagnet inter- 
faces are characterized by the dimensionless parameters Jbx 
and 7b2, respectively, where 7^1.^2 = RBi,B2<yn/^n, Rbi.b2 are 
the resistances of left and right S/F interfaces, respectively, 
(7„ is the conductivity of the F layer, ^„ = ^JdJ/T/kTc, Df 
is the diffusion coefficient in the ferromagnetic metal and 7^ 
is the critical temperature of the superconductor (we assume 
h = kB = 1). We also assume that the S/F interfaces are not 
magnetically active. We will consider diffusive limit, in which 
the elastic scattering length I is much smaller than the decay 
characteristic length ^fi. In this paper we concentrate on the 
case of a SIFS tunnel Josephson junction, when 7^1 ^ 1 (tun- 
nel barrier) and 7^2 = (fully transparent interface). For com- 
parison, we also consider two other limiting cases: an SFS 
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-d/2 df/2 



FIG. 1: Geometry of the considered system. The thickness of the 
ferromagnetic interlayer is df. The transparency of the left S/F inter- 
face is characterized by the Jbx coefficient, and the transparency of 
the right F/S interface is characterized by y^l- 



junction (7bi = Jbi = 0) and a SIFIS junction (Ybi,7b2 > 1). 

Under conditions described above, the calculation of the 

Josephson current requires solution of the one-dimensional 

Usadel equations, In the F layer the equations has the 
forrr22=2i 



D 



[{(0 ± ih)a, + X— o-zG/Ta)(^z, 



(1) 



where positive sign ahead of h corresponds to the spin up state 
(t) and negative sign to the spin down state (|), (O — 2nT{n + 

are the Matsubara frequencies, h is the exchange field in the 
ferromagnet and a, is the Pauli matrix in the Nambu space. 
The parameter T,„ is the spin-flip scattering time. The influ- 
ence of spin-flip scattering on various properties of S/F struc- 
tures was considered in a number of papers . ^ ^t^^?- ^ 1^2,23^3 1^32 
We consider the ferromagnet with strong uniaxial anisotropy, 
in which case the magnetic scattering does not couple the spin 
up and spin down electron populations. 

The Usadel equation in the S layer can be written a&i^ 



(2) 



where Ds is the diffusion coefficient in the superconductor. In 
Eq. (|2|l Gs = and we omit subscripts 'f (j)' because 

equations in superconductor look identically for spin up and 
spin down electron states. 

In Eqs. ([T]i-(|2|i we use following matrix notations (we omit 
'f, 's' and 't (i)' subscripts) 



G(jc,a)) 



G 

p* 



A(x) = 



^{x) 

A*(x) 



(3) 



where G and F are normal and anomalous Green's functions, 
respectively, and A(x) is the superconducting pair potential. 
The matrix Green's function G satisfies the normalization con- 
dition. 



1, G-+FF* = \, 



(4) 



and the pair potential ^(x) is determined by the self- 
consistency equation 



(i)>0 



(5) 



The boundary conditions for the Usadel equations at the left 
and right sides of each S/F interface are given by relations^'* 



^ny[G,-G, 



2<^„7m ( Gf—Gf 



±dfl2 



-dfll 



±dfl2 



2^„7.2 [Gf^Gfj^^^=[Gf,Gl,^i^, 



(6a) 
(6b) 
(6c) 



where y ~ £,sG„/£,„Gs, C7s is the conductivity of the S layer and 

To complete the boundary problem we also set boundary 
conditions at x — ±oo. 



G,(±-) = 



0) 



v/|A|2 + «2' 



|Ak"P/2 



(0^ 



CO-" 



(7a) 



F, (7b) 



where (p is the superconducting phase difference between S 
electrodes. 

In Matsubara technique it is convenient to parameterize the 
Green's function in the following way, making use of the nor- 
malization condition, Eq. 



cos sin Oe''^ 
sinOe^'^ — cos0 



(8) 



Solving a system of nonlinear differential equations, 
Eqs. ([T]i-0^ generally can be fulfilled only numerically. We 
present full numerical calculation in Sec. |IV] The analytical 
solution can be constructed in case of one S/F bilayer, when 
we can set the phase X in Eq. ([8]) to zero. We can also set 
the phase to zero in case of long S/F/S junction, where the 
thickness of the feiTomagnetic layer df ^ ^fi. In that case, 
the decay of the Cooper pair wave function in first approxi- 
mation occurs independently near each interface. Therefore 
we can consider the behavior of the anomalous Green's func- 
tion near each S/F interface, assuming that the ferromagnetic 
interlayer is infinite. This analytical calculation for an S/F/S 
trilayer with long ferromagnetic interlayer is performed in the 
next section. 

The general expression for the supercurrent is given by 
i7tT(7n +^ f- d d ~ \ 

E ^\'f<^d-/f^-^f<^Tx'f<^)^ (9) 

where F^m)(x,a)) = Fj^^^^{x,—(X)) are the anomalous Green's 
functions in the ferromagnet. 
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m. CRITICAL CURRENT OF JUNCTIONS WITH LONG 
FERROMAGNETIC INTERLAYER 

We need to solve the complete nonlinear Usadel equations 
in the ferromagnet, Eqs. ([U. For SIFS junctions, an analytical 
solution may be found if df ^ ^/i and we can set the phase 
of the anomalous Green's function to zero (see discussion in 
SeclID. 

Setting Xs = Xf = we have the following 9- 
parameterizations of the normal and anomalous Green's func- 
tions, Eq. ([8]), G = cos 9 and F ~ sin9. In this case we can 
write Eqs. ^ in the F layer as 



(0±ih- 



cos 9 



'nil) 



sin 9 



nil)- 



(10) 



2 dx^ \ T„] 
In the S layer the Usadel equation, Eq. may be now written 



where £,/ = \jDfjh and the boundary condition 9f(x^^)- 
has been used. In Eq. ( fT6l ) we use the following notations 



q = y/ljhyjoi ± ;•/! + 1 /t,„, 
= (1/t,„) (a)± ;•/! + ! /t,„)" 



(17a) 
(17b) 



Here we again adopt convention that positive sign ahead of 
h corresponds to the spin up state (t) and negative sign to 
the spin down state (|). Here and below we did not write 
spin labels 't (i)' explicitly but imply them everywhere they 
needed. 

For the right interface (x = df/2), a first integral of Eq. (fTOl ) 
leads to a similar equation. 



2 dx 



• % /i 9-2 % 

asm — W 1 — e^^sm — . 
^ 2 V 2 



(18) 



Following Faure et alM we integrate Eq. ( fT6T l, which gives 



Ds d^9s 



(0 sin 9s — A(x) cos 9s 



(11) 



2 dx^ 

The self-consistency equation in the S layer acquires the form 
A(x)ln^ = 7rr^ ("^^-sin^-sine^V (12) 



(U>0 



In the case of Xs = = 0, the boundary conditions, 
Eqs. (|6]l, for the functions 9f ,s at each S/F interface can be 
written as 



V '^■^^ / ±df/2 / ±dfll 



39 



^'.Ybi(^] = sin {9f- 9s) 



-df/2 



39 

^„YB2[^)^^^^^sm{9s~9f] 



-df/2 ' 



df/2 ■ 



(13a) 
(13b) 
(13c) 



■ e^sm -j- — cos -j- 



1 — g'^sm -4- +COS -d- 



= giexp -2q 



df/2- 



(19) 



The integration constant gi in Eq. (fT9] l should be deter- 
mined from the boundary condition at the left S/F interface, 
Eq. ( |13bt . Since we consider the tunnel limit (7^1 ^ 1), we 
can neglect small 9f in the right hand side of Eq. ( I13bb and 
also assume, neglecting the inverse proximity effect. 



9s{—df/2) = arctan — . 



Then Eq. ( I13bb becomes 



d9f 



-df/2 



'G{n), G(n) : 



(20) 



(21) 



From Eqs. (fTST i and (ISTT i we obtain the boundary value of 9f 
six = —df/2 and substituting it into Eq. (fT9l ) we finally get 



The boundary conditions at x = are 



a (±00) ^ arctan — 



(14) 



In the equation for the supercurrent, Eq. (|9l), the summation 
goes over all Matsubara frequencies. It is possible to rewrite 
the sum only over positive Matsubara frequencies due to the 
symmetry relation. 



%)T(«) = %(.v)i(-«)- 



(15) 



In what follows, we will use only o > in equations contain- 
ing (O. 

For the left interface (tunnel barrier at x = —df/2), a first 
integral of Eq. ( fTOl i leads to 



^1 



G^in) 1 



I67J1 



(22) 



Linearizing Eq. ( fT9] l, we can now obtain the anomalous 
Green's function in the ferromagnetic layer of the SIF tun- 
nel junction with infinite F layer thickness. Similar formula 
for the FS bilayer with a transparent interface (7^2 = 0) was 
developed by Faure et al?^ [to obtain it one should inte- 
grate Eq. ( fTSl ) and then linearize the resulting equation]. The 
anomalous Green's function at the center of the F layer in a 
SIFS junction may be taken as the superposition of the two 
decaying functions, taking into account the phase difference 
in each superconducting electrode. 



9f = 



^2exp q 



/^lexp -q 



'df/2 ,(p 



df/2- 



-X (p 
--'2 



(23) 
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FIG. 2: (Color online) The F layer thickness dependence of the 
critical current for SFS (7^1,2 = 0), SIFS (7^1 = 10^, Ybi = 0) and 
SIFIS (7^1.2 = 10^) junctions in the absence of spin-flip scattering. 
Red dashed lines correspond to the modulus of the analytical results 
J31b.(l25t and l |29l l and black solid lines correspond to the results of 
numerical calculation in Sec.|lV] h = 37lTc, T = 0.5Tc. 



The expression for g2 was obtained in Ref. |23| for the 
rigid boundary conditions at the transparent FS interface, 
9f{df/2) = arctan (|A|/a)) and reads 



g2 



F{n) = 



(l-£2)F2(n) 



W{l-e^)F\n)- 
|A| 

(£)+ y/W+\W' 



1 + 1]2 



(24a) 



(24b) 



Using the above solutions and Eqs. ( fTSl ) we arrive at sinu- 
soidal current-phase relation in a SIFS tunnel Josephson junc- 
tion with the critical current 



16kT , 



■Re 



G{n)F{n)Qxp{-qdf/t,f) 
1 ^J{l-e^)F\n) + \ + l 



(25) 



definition 

2 _ 



Here and below we fix positive sign in the 
of q, in Eqs. ( fTTI ): q — ^Ijh^ Oi + ih + 1/t, 
(1/ T„,) (o + ih + 1/ T„,)^' . It is possible since we already per- 
formed summation over spin states and have to define now 
spin-independent values. In Eq. (l25l l and below Rfi is a full 
resistance of an S/F/S trilayer, which include both interface 
resistances of left and right interfaces and the resistance of the 
ferromagnetic interlayer. In case of SIFS and SIFIS junctions 
the F layer resistance can be neglected compared to large re- 
sistance of the tunnel barrier. 

At this point we define the characteristic lengths of the de- 
cay and oscillations ^f\2 as, 



q/^f=\l^fX+i/^f2, 



(26a) 



FIG. 3: (Color online) The F layer thickness dependence of the criti- 
cal current in a SIFS junction [modulus of the Eq. l l25b 1 for different 
values of a = 1 jnTctm, h = 37tTc, T = O.STc. 



The critical current in Eq. (I25b is proportional to the small 
exponent s\p(^—df/^fi). The terms neglected in our ap- 
proach are of the order of exp (— 2ii^/^^i) and they give a 
tiny second-harmonic term in the current-phase relation. 

The critical current equation ( l25T l can be simplified in the 
limit of vanishing magnetic scattering, T,^ ' <C TzTc, 



11=0 



G(n)F(n)exp (^-^^ cos 
Vf2(«) + 1 + 1 



(27) 



Eq. ( l25T l also simplifies near Tc and may be written as (for 

Tc < h) 



TtA^ { df 



df 

172 



■ (28) 



The damped oscillatory behavior of the critical current can be 
clearly seen from this equation. With increasing df the junc- 
tion undergoes the sequence of Q-n transitions when positive 
values of the IcRn product correspond to a zero state and neg- 
ative values coiTespond to a ;r state. 

Eq. ( |28] l in the absence of spin-flip scattering coincides with 
the corresponding equation, Eq. (37), from the Ref. 17, taken 
in the limit of long df ^ ^fi in case of 7gi » 1, 752 = 0. 

Using the same approach we can obtain the equation for 
the critical current in a SIFIS structure with two strong tunnel 
baiTiers between the feiTomagnet and both superconducting 
layers (7^1,2 > 1), 



eq„ Yb\7b2 



G2(„)exp(^) 



q 



(29) 



1 



This formula coincides with corresponding expression 
Eq. (39) for the critical current in a SIFIS structure in Ref.l23l 
for 7bi^2 = 7b >• 1 and df ^ t,f\. Eq. ( |29] l near Tc may be 
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FIG. 4: (Color online) The F layer thickness dependence of the crit- 
ical current in a SIFS junction [modulus of the Eq. J25H for differ- 
ent values of exchange field h in the absence of spin-flip scattering, 
T = Q.5Tr. 



FIG. 5: (Color online) Fit to the experimental data from Ref. 14 for 
the critical current in a Nb/Al203/Nio.6Cuo.4/Nb junction. The fitting 
parameters are: h/kg = 950 K and 1 /%„ = 1.6 h. 



written as (for ^ h) 

n\Af-£,f2 Ybi + Ybi 



IcRn — 



2eTc^„ YbiYbi 



X cos (^)exp( Isinl^-^^ 



in 



if2 



(30) 



where is defined by tan(^) ~ (^^2/^/1 ■ Eq. dSOl l in the ab- 
sence of spin-flip scattering coincides with the corresponding 
equation, Eq. (35), from the Ref.[T3, taken in the limit of long 

We also provide here equation for the critical current in an 
SFS junction [see Ref. [23, Eq. (74)], written in our notations. 



IcRn 



6AnTd 



Re 



~ F^{n)qtx^{-qdf/E,t) 
hW {I -£^)F^{n) + 1 + 1? 



(31) 



We compare critical current dependencies over dj for SFS 
[Eq. (EBl, SIFS [Eq. (|25]l] and SIFIS [Eq. ^] structures in 
Fig.|2] Each of above junction types undergoes the sequence 
of Q-n transitions with increasing thickness of the F layer 
From the figure we see that the transition from to ;r state 
occurs in SIFS tunnel junctions at shorter df than in SFS junc- 
tions with transparent interfaces, but at longer df than in SIFIS 
junctions with two strong tunnel barriers. This tendency can 
be qualitatively explained by the fact that in structures with 
barriers (SIFS, SIFIS) part of the n phase shift occurs across 
the barriers. Therefore a thinner F layer in a SIFS junction 
compared to an SFS one is needed to provide the total shift of 
n due to the order parameter oscillation. For the same reason, 
O-TT transition in a SIFIS junction occurs at a smaller thick- 
ness than in a SIFS junction. We note that in Fig. |2] we plot 
both analytical and numerical calculated Ic{df) dependencies, 
where numerical calculation was performed for full boundary 
problem Eqs. ([T]i-(|7]i [see further discussion in Sec. lIVII . 



In Fig. [3] we plot the F layer thickness dependence of the 
critical current in a SIFS junction for different values of spin- 
flip scattering time. For stronger spin-flip scattering the pe- 
riod of supercurrent oscillations increases and the point of Q-n 
transition shifts to the region of larger df. The same tendency 
exists for SFS and SIFIS junctions. 

In Fig. m we plot the F layer thickness dependence of 
the critical current in a SIFS junction for different values of 
the exchange field h. We see that for large exchange fields 
h ^ nTc the critical current scales with the ferromagnetic co- 
herence length ^f. 

From comparison with numerical results presented in Fig.|2] 
we can conclude that the results for the critical current in SIFS 
junctions presented in Figs.[3]|4]give correct magnitude of the 
IcRn product for df > £,„/2. 

As an application of the developed formalism, we present 
in Fig. |5] the theoretical fit of the experimental data for a 
Nb/Al203/Nio.6Cuo.4/Nb junctions by Weides et al}* mak- 
ing use of Eq. ( l25T l. We used following values of parame- 
ters: Rb = 3.9 m£2, Df = 3.9 cm^/s, T = 4.2 K}'^ Tc = 7.2 K 
(damped critical temperature in Nb). Good agreement was ob- 
tained with the following parameters: h/kB = 950 K, l/tm = 
1.6 h (see Fig.|5]l. These parameters can be compared with 
parameters obtained by Oboznov et alM for similar ferro- 
magnetic material, Nio.53Cuo.47: h/kB = 850 /T, 1 /%„ — 1.3 h. 
Higher Ni concentration in the NiCu alloy in the experiment 
of Weides et al. results in higher exchange field. 

In Ref. it was suggested that a "dead" layer exists in the 
ferromagnet near each S/F interface, which does not take part 
in the "oscillating" superconductivity. Other authors also in- 
clude into consideration the existence of nonmagnetic layers 
at the interface of the ferromagnet and the superconductor or 
normal meta L^^'^^i^^ Thickness of the "dead" layer cannot be 
calculated quantitatively in the framework of our model and 
also can not be directly estimated from the experiment. In 
the experiment of Weides et al. the range of F layer thick- 
nesses was rather narrow and only the first Q-n transition was 
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FIG. 6: DOS on the free boundary of the F layer in the FS bilayer 
calculated numerically in the absence of spin-flip scattering for dif- 
ferent values of the F layer thickness df. (E) (dashed line), (E) 
(dotted line) and N(E) (solid line), £„■ = SkTc, T = O.STc. (a): 
df/^„ = 0.4, (b): df/^„ = 1, (c): df/^„ = 1.6, (d): df/^„ = 2.2. 



FIG. 7: DOS N{E) on the free boundary of the F layer in the FS 
bilayer calculated numerically for a = 1/ TzTc'Cm = (solid line), a = 
0.5 (dashed line) and a = I (dotted line) for different values of the 
F layer thickness df, Eex = ^itTc, T = Q.STc. (a): df/t,„ = 0.4, (b): 
df/t,,, = 1, (c): df/^n = 1.6, (d): df/t,n = 2.2. 



observed. Due to these reasons we did not take into account 
the existence of a nonmagnetic layer in our fit. This question 
deserves separate detail experimental and theoretical study. 

We should mention that the above estimates of exchange 
field and spin-flip scattering time could be different if we 
consider magnetically active S/F interfaces. It was shown in 
Ref . |28 that the effect of spin-dependent boundary conditions 
on the superconducting proximity effect in a diffusive ferro- 
magnet results in the change of the period of critical current 
oscillations. 



rv. CRITICAL CURRENT OF JUNCTIONS WITH 
ARBITRARY LENGTH OF THE FERROMAGNETIC 
INTERLAYER 

In the previous section we derived the expression for the 
critical current of a SIFS junction in case of considerably long 
F layer thickness, df^£,fi. For arbitrary F layer thickness in 
the absence of spin-flip scattering, general boundary problem 
([ril-(|7]l was solved numerically using the iterative procedure. 
Starting from trial values of the complex pair potential A(x) 
and the Green's functions Gs.f, we solve the resulting bound- 
ary problem. After this we recalculate Gsj and A(x). We 
repeat the iterations until convergency is reached. The self- 
consistency of calculations is checked by the condition of con- 
servation of the supercurrent across the junction. 

In Fig. |2] we compare numerically and analytically calcu- 



lated Ic{df) dependencies in case of SFS, SIFS and SIFIS 
junctions. We see that, as expected, the numerical method pro- 
vides correction only for small length of ferromagnetic layer. 
We note that for SFS and SIFS junctions analytical curves dSTT l 
and (|25] | practically coincide with numerical results in the re- 
gion of the first O-TT transition. For a SIFIS junction this transi- 
tion occurs at smaller df, where the assumptions of the section 
Uniare not valid. However, in presence of strong spin-flip scat- 
tering the first O-TT transition peak in a SIFIS junction shifts to 
the region of larger df and Eq. ( |29] | describes the transition 
accurately. 

The main result of this section is that Eq. dZSl ) for the crit- 
ical current of a SIFS junction can be used as a tool to fit 
experimental data in SIFS junctions with good accuracy. 



V. DENSITY OF STATES OSCILLATIONS IN THE 
FERROMAGNETIC INTERLAYER 

It is known that in a ferromagnetic metal attached to the 
superconductor the quasiparticle DOS at energies close to the 
Fermi energy has a damped oscillatory behavior^^'^"*'^^ Ex- 
perimental evidence for such behavior was provided by Kon- 
tos et ali^ In SIFS junctions we can compare the DOS oscil- 
lations with the critical current oscillations. 

We are interested in the quasiparticle DOS in the F layer 
in the vicinity of the tunnel barrier (x — —df/2 + in Fig. [T]). 
Below we will refer to the local DOS at this point. For the case 
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FIG. 8: (Color online) The F-layer dependence of the function 
SN{df) in the absence of spin-flip scattering, h = 37tTc , T = 0.5Tc. 
Black solid line is a result of the numerical calculation; blue dashed 
line is calculated with the use of Eq. ( 141b . Red line shows normalized 
critical current for a SIPS junction. Zero and K states defined from Ic 
are indicated by red color, while zero and n states defined from the 
DOS are indicated by black color. 



of strong tunnel barrier (7^1 ^ 1) left S layer and right FS bi- 
layer in Fig. [T] are uncoupled. Therefore we need to calculate 
the DOS in the FS bilayer at the free boundary of the fer- 
romagnet. Solving numerically Eqs. (fT0ll- (fl4l l. we set to zero 
the Of derivative at the free edge of the FS bilayer, x = —df/2, 

We use the self-consistent two step iterative 
procedure^Si^Sili. In the first step we calculate the pair 
potential coordinate dependence A(x) using the self- 
consistency equation in the S layer, Eq. ( fT2] i. Then, by 
proceeding to the analytical continuation in Eqs. (fTOl i. (fTTT i 
over the quasiparticle energy ico ^ E + iO and using the A(jic) 
dependence obtained in the previous step, we find the Green's 
functions by repeating the iterations until convergency is 
reached. We define the full DOS N{E) and the spin resolved 
DOS A^m) (£■), normalized to the DOS in the normal state, as 

NiE)=[N^iE)+Ni{E)]/2, (32a) 
N^i) (£) = Re [cos 9^ ^) {i(0 E + iO)] . (32b) 

The numerically obtained energy dependencies of the DOS 
at the free F boundary of the FS bilayer are presented in Figs.|6] 
and|7] Fig. |6] demonstrates the DOS energy dependence for 
different df.Ai small df full DOS turns to zero inside a mini- 



gap, which vanishes with the increase of df. Then the DOS at 
the Fermi energy A^(0) rapidly increases to the values larger 
than unity and with further increase of df it oscillates around 
unity while it's absolute value exponentially approaches unity 
(see also Fig.[8]l. In Fig.|6]we also plot the spin resolved DOS 
energy dependencies Aff(£') andA^|(£'). Fig. |7] demonstrates 
full DOS energy dependence for different values of spin-flip 
scattering time. For stronger spin-flip scattering the minigap 
closes at smaller df, the period of the DOS oscillations at the 



Fermi energy increases and the damped exponential decay oc- 
curs faster. 

In case of long F layer {df 3> <^/i) it is also possible to ob- 
tain an analytical expression for the DOS at the free boundary 
of the ferromagnet, 



1 2 

1--Re0^^^), 



(33) 



where Ob]{[) is ^ boundary value of Qf at x = —df/2. It can 
be obtained by the mapping method, similar to the one used in 
the electrostatic problems. We consider the FS bilayer where 
X G [—df/2, df/2] stands for the ferromagnetic metal and 
x> df/2 stands for the superconductor; the point x = —df/2 
corresponds to the free F layer boundary. For infinite F layer 
{df °°) the solution for Qf^(i) far from the interface is given 
by the exponential term in Eq. ( |23] ). written in the real energy 
space. 



e 



fUi) 



4 ^ ( x-df/2 \ 



where 



82 = 
FiE) = 



y^2/hy^-iER ± ih + 1 / T,„, 

il/T„,)i-iER±ih+l/T,„)-\ 

{l-n^)F\E) 

[V(1-J?2)F2(£) + 1 + 1]2' 

|A| 



Er=E + iO. 



(34) 

(35a) 
(35b) 

(35c) 
(35d) 



-iER+J\A\^-Ej 



Here, as above, positive sign ahead of h corresponds to the 
spin up state in Eq. (|34|) and negative sign for the spin down 
state. By using the arrow 'from right to left' in 9 we 
want to stress that this solution is induced in the ferromagnet 
from the right FS interface. 

In the case of finite ferromagnet length the boundary con- 
ditions at the free F layer boundary, x = —df/2, become 



de 



/T(i) 



: 0. (36) 



-df/2 



To ensure these conditions we add another exponential solu- 
tion. 



'/TU) 



^2exp 



3df/2- 



(37) 



resulting from the mirror image of the F layer with respect 
to the point x = —df/2. At x — —df/2 both exponential 
terms are equal to each other and the final solution, = 

9 1///2)+ 8 (—t/^/2), is two times largerthan the 

solution for infinite ferromagnetic layer at this point and reads 



9, 



8F{E) 



V(1-772)F2(£) + 1 + 1 



exp 



17. 



(38) 



8 



This equation coincides with the resuh obtained in Ref . 112 by 
direct integration of the Usadel equation. 

In Fig. m we plot analytically and numerically calculated 
function 



5N{df) = \l-No\, No=N{E = 0), 



(39) 



together with the Ic{df) dependence for a SIFS junction. We 
see that the point of Q-n transition on the Ic{df) plot does not 
coincide with the first minimum of 5N{df) corresponding to 
sign change of 1 — A^o- This difference can be qualitatively ex- 
plained as follows. The transition from to tt state in a junc- 
tion, seen as sign change of Ic{df), is the result of interference 
of solutions for Qf originating from two S electrodes. Q-n 
transition in Ic{df) occurs approximately at such thickness df 
when the boundary value of Qf in Eq. ( |23T l at jc = —df/2 be- 
comes negative, i.e. when df acquires the phase shift n. On 
the other hand, sign change of 1 — A^o occurs at such df when 
the boundary value Qi, in Eq. ( |38] | becomes an imaginary num- 
ber, i.e. when Qf acquires the phase shift n/2. It occures at 
smaller df compared to O-tt transition in the critical current. 
Corresponding and K states defined from Ic and from the 
DOS are indicated in Fig. |8] 

It is also seen from Fig.[8]that the DOS oscillations have the 
period approximately twice smaller than those of the critical 
current. This fact is easy to see from the analytical expression 
for dN{df). Using Eqs. (O we obtain 



5N{df) = 32 



Re 



1 



(./2- + 1)2 



exp -po- 



2df 

17 



(40) 



where tjq = TJ (£ = 0) and po — p{E~ 
At vanishing magnetic scattering, T„ 
can be simplified. 



= 0) in Eqs. (I35ali-(l35b]). 
^ ^ TiTc, this equation 



VI. CONCLUSION 



We have developed a quantitative model, which describes 
the oscillations of the critical current as a function of the F 
layer thickness in a SIFS tunnel junctions with thick ferro- 
magnetic interlayer, df » ^/i, in the dirty Hmit. We justified 
this model by numerical calculations in general case of arbi- 
trary df: for all values of parameters characterizing material 
properties of the ferromagnetic metal numerical and analyti- 
cal results coincide in physically important region of the first 
Q-7Z transition. Thus the derived analytical expression for the 
critical current can be used as a tool to fit experimental data in 
various types of SIFS junctions. We have discussed the details 
of the damped oscillatory behavior of the critical current for 
different values of the F layer parameters. 

We also studied the superconducting DOS induced in a fer- 
romagnet by the proximity effect. We showed that the oscilla- 
tion pattern of DOS at the Fermi energy in the ferromagnet (at 
location of the tunnel junction) does not coincide with that of 
the critical current in a SIFS junction and it's period is approx- 
imately twice smaller Therefore the DOS oscillations do not 
reflect the Q-n transition in Ic{df). We calculated the quasi- 
particle DOS in the F layer in the close vicinity of the tunnel 
barrier which can be used to obtain current-voltage character- 
istics for a SIFS junction. These calculations wiU be presented 
elsewhere. 

Finally, we used our results to fit recent experimental data 
for SIFS tunnel junctions and extracted important parameters 
of the ferromagnetic interlayer. 



5N{df) 



32 



3+2V2 



exp 



2df\ /2df 



(41) 



where characteristic lengths of decay and oscillations ^fi2 are 
given by Eq. ( |26b| i with the substitution i(0 ^ E + iO. This 
equation can be compared with Eq. ( |27] |. We see that the pe- 
riod of the DOS oscillations is approximately twice smaller 
than the period of the critical current oscillations and the ex- 
ponential decay is approximately twice faster than the decay 
of the critical current. 
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